Three-body Casimir effects and non-monotonic forces 
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Casimir interactions are not pair-wise additive. This property leads to collective effects that we 
study for a pair of objects near a conducting wall. We employ a scattering approach to compute 
the interaction in terms of fluctuating multipoles. The wall can lead to a non-monotonic force 
between the objects. For two atoms with anisotropic electric and magnetic dipole polarizabilities 
we demonstrate that this non-monotonic effect results from a competition between two- and three 
body interactions. By including higher order multipoles we obtain the force between two macroscopic 
metallic spheres for a wide range of sphere separations and distances to the wall. 

PACS numbers: 42.25.Fx, 03.70.+k, 12.20.-m. 



A hallmark property of dispersion forces is their non- 
additivity which clearly distinguishes them from elec- 
tromagnetic forces between charged particles [T]. Work 
on the interactions between multiple objects is limited 
mostly to atoms or small particles which are described 
well in dipole approximation [2]. This approximation 
cannot be used for macroscopic objects at separations 
that are comparable to their size since higher order 
multipole fluctuations have to be included |31 Sj- In 
such situations, also other common "additive" methods 
such as proximity or two-body-interaction approxima- 
tions fail. Three-body effects for macroscopic bodies have 
been studied in quasi two-dimensional (2D) geometries 
that are composed of parallel perfect metal cylinders of 
quadratic [S] or circular [HI [7] cross section and parallel 
sidewalls. For this setup non-monotonic forces have been 
found and interpreted as resulting from a competition 
between electric and magnetic polarizations which are 
decoupled for quasi 2D geometries of perfect metal struc- 
tures. In this Letter we investigate collective 3-body ef- 
fects between compact objects, including anisotropic po- 
larizabilities, and a wall in three dimensions using a re- 
cently developed scattering approach |11[5]. This allows 
us to observe the influence of polarization coupling and 
anisotropy on non-monotonic effects. 

We consider the retarded Casimir interaction between 
a pair of atoms with anisotropic electric and magnetic 
polarizabilities near a conducting wall, see Fig. [l] We 
identify a competition between 2- and 3-body effects and 
prove that this leads to a non-monotonic dependence 
of the force between the atoms on the wall separation 
H for each of the four possible polarizations of fluctu- 
ations (electric/magnetic and parallel/perpendicular to 
the wall) separately. For isotropic polarizabilities we find 
that only the force component due to electric fluctuations 
is non-monotonic in H . This flndings suggest the possi- 



bility to engineer the monotonicity properties of the force 
by suitable tailoring of the polarization tensors via the 
shape and material composition of macroscopic particles. 

For atoms, magnetic effects are almost always rather 
small in the retarded limit. Contrary to this, for con- 
ducting macroscopic objects contributions from electric 
and magnetic multipole fluctuations are comparable. To 
study the effect of higher-order multipoles, we consider 
also two perfect metal spheres near a wall, see Fig. [T| 
Based on consistent analytical results for large separa- 
tions and numerical computations at smaller distances 
we find a non-monotonic dependence of the force between 
the spheres on H. Unlike for atoms, this effect is occurs 
at sufficiently large sphere separations only. 




FIG. 1: Geometry of the two-sphere/atom and sidewall sys- 
tem. Shown are also the mirror images (grey) and two- and 
three-body contributions (solid and dashed curly lines, re- 
spectively) . 

As derived by Refs. |4t jSj, the Casimir energy of two 
bodies in the presence of a perfectly conducting sidewall 
can be obtained within a scattering approach by employ- 
ing the method of images that introduces ffuctuating cur- 
rents on the mirror bodies. The Casimir energy of the 
original system is then given by the energy of the origi- 
nal and the image objects and it can be expressed as an 
integral over imaginary wave number, 
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which is given by the T-matrix T that relates the incom- 
ming and scattered electromagnetic (EM) fields for each 
body, and by the U-matrices U"'', U^'"'' that describe 
the interaction between the multipoles of object a and 
object (3 and between the multipoles of object a and the 
image of object f3, respectively. The T-matrix depends 
only on properties of the bodies such as polarizability 
or size and shape. The U-matrices depend only on the 
distance vector between the objects and decay exponen- 
tially with distance and wave number k. The matrix Moo 
accounts for the subtraction of the object's self-energies 
and hence follows from M by taking the limit of infinite 
separations, i.e., by setting all U-matrices to zero. For a 
multipole expansion the matrix elements are computed in 
a vector spherical basis for the EM field with partial wave 



numbers I > I, m = —I, . . . ,1. (Details of this expansion 
and the U-matrix elements can be found in Ref. [8j.) 
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FIG. 2: Typical orientations of electric (E) and magnetic (M) 
dipoles and image dipoles for H/L — > and H/L — > oo. 

In the following we study the force F — —d£/dL be- 
tween the two objects at separation L and hence elim- 
inate the contributions to the energy that depend only 
on the sidewall separation H, see Fig.[T] We expand the 
determinant of Eq. as 



det(MM;^i) = det(l + TU^) det(l + TU^) det [l - (1 + TV')-^T{V^^ + V^'^^){1 + TV^)-'^T{V^^ + V^'^^)] . (3) 

The first two determinants on the r.h.s. yield together twice the interaction energy between a single object and the 
sidewall since = U^'^^ = U^'^^ describes the multipole coupling between one object and its image and hence 
depends only on H. Hence, we consider only the energy £oo that corresponds to the last determinant of Eq. ([s]) and 
provides the potential energy of the two objects in the presence of the sidewall so that F = —d£oo/dL. In the absence 
of the sidewall, iJ — > oo, the matrices U^'"'' all vanish and £oo simplifies to the energy between two spheres [4 . For 
an interpretation in terms of multiple scatterings, it is instructive to use the relation In det = Trln and to expand the 
logarithm which yields 
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where we have written the inverse matrices of Eq. ([3| 
as series. The trace acts on an alternating product of T- 
and U-matrices which describe scattering and free propa- 
gation of EM fluctuations, respectively. Multiple scatter- 
ings between an object and its image (TU^) are followed 
by a propagation to the other object-image pair, either 
to the object (U^^) or its image (U^'^i)^ between which 
again multiple scatterings occur before the fluctuations 
are scattered back to the initial object or its image (U^^ 
or and the process repeats. This expansion is use- 

ful for small objects or large separations. 

As the first application, we consider the case of two 
identical, ground state atoms near a wall, see Fig. [T| The 
separation between the atoms is L and they have equal 
distance H from the wall. In dipole approximation, the 
retarded limit of the interaction is described by the static 
electric (a^, ay) and magnetic {Pz, dipole polarizabil- 
ities of the atoms which can be different perpendicular (z) 
and parallel (||) to the wall. The T-matrix of the atoms 
is diagonal and has finite elements only for the dipole 
channel (partial waves with ^ = 1), given by T^^ = ^azK^, 



^ for electric and 



for magnetic polarization with m — ±1. For atoms, the 
polarizability is much smaller than L^, and hence it is 
sufficient to compute the interaction to second order in 
the polarizabilities. This amounts to consider only the 
terms with p — 1, n = n' = Oof Eq. The resulting 
energy f oo is then compared to the well-known Casimir- 
Polder (CP) interaction between two atoms (without the 
wall) 
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which corresponds to the sequence TU^-'^TILJ-'^^ jn Eq. Q. 
The total interaction energy is 

£oo_{L, H) = S2,\{L) + 52a(^. L) + S^{D, L) (6) 

with D = v^JX+UP. The 2-body energy £2^\{D, L) 
comes from the sequence TU^'^^TU^'^^ in Eq, Q and 
hence is the usual CP interaction between one atom and 
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the image of the other atom (see Fig. [T]). The change in 
the relative orientation of the atoms with t = L/D leads 



to the modified CP potential 
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The 3-body energy £3(0, L) corresponds to the sequences TU^iTU^'i^ TU^'2iTUi2 in Eq. Q and hence describes 
the collective interaction between the two atoms and one image atom. It is given by 
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For isotropic electric polarizable atoms this result agrees 
with that of Ref . |2] . It is instructive to consider the two 
hmits H <^ L and H :^ L. For H L one has D — > 
L and the 2-body potentials are identical, £2.\{L,L) — 
£2,1 (i)- The 3-body energy becomes 



£s{L,L) 



he 



-33a| -I- 13a^ + 14q!||/3^ - {a ^ (i) 

(9) 

The total energy f 00 is now twice the energy of Eq. ([s]) 
plus the energy of Eq. Q and hence foo becomes the 
CP potential of Eq. ([s]) with the replacements Uz — *■ , 
a\\ ^ 0, Pz ^ 0, /3|| 2/3||, i.e., the 2-body and 3- 
body contributions add constructively or destructively, 
depending on the relative orientation of a dipole and its 
image which together form a dipole of zero or twice the 
original strength (see Fig. |2]). For H ^ L the leading 
correction to the CP potential of Eq. Q comes from the 
3-body energy which in this limit becomes to order H^^ 
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The signs of the polarizabilities in the leading term 
~ H~'^ can be understood from the relative orientation of 
the dipole of one atom and the image dipole of the other 
atom, see Fig. [2] If these two electric (magnetic) dipoles 
are almost perpendicular to their distance vector they 
contribute attractively (repulsively) to the potential be- 
tween the two original atoms. If these electric (magnetic) 
dipoles are almost parallel to their distance vector they 
yield a repulsive (attractive) contribution. For isotropic 



polarizabilities the leading term of Eq. ( 10 ) vanishes and 
the electric (magnetic) part ~ of the 3-body energy 
is always repulsive (attractive). 

The above results determine the variation of the force 
between the two particles with H . If the two particles 
have only either or /?|| polarizabihty, their attractive 
force is redueed when they approach the wall from large 
H due to the repulsive 3-body energy. At close proximity 
to the wall the fluctuations of the dipole and its image 
add up to yield a force between the particles that is en- 
haneed by a factor of 4 compared to the force for H ^ 00. 
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Equivalent arguments show that the force between par- 
ticles with either ay or f3z polarizability is enhaneed at 
large H and redueed to zero for H 0. This proves 
that the force between particles which both have either 
of the four polarizabilities is always non-monotonic. The 
situation can be different if more than one polarizability 
is finite, especially for isotropic particles. In the latter 
case all contributions (electric, magnetic, mixed) are en- 
hanced for iJ — > and only the electric term is reduced 
at large H so that only the electric part gives a non- 
monotonic force. In general, the monotonicity property 
depends on the relative strength and anisotropy of the 
electric and magnetic polarizabilities. 

As the second application we study two macroscopic 
perfect metallic spheres of radius R for the same 
geometry as in the case of atoms where the lengths 
L and H are measured now from the centers of the 
spheres, see Fig. [T] The T-matrix is diagonal and the 
elements = (-l)'f /z+i/2(«i?)/i^z+i/2(Ati?), T^^ = 

(-l)'f[/, + l/2(A«i?,) + 2/ti?/;^^/2(«i?)]/[AVl/2(«i?) + 

2kRKi^^^^{kR)] are given in terms of the Bessel func- 
tions I^, K^. First, we expand the energy in powers 
of R by using Eq. ^ which implies that we expand 
the T-matrices for small frequencies but use the exact 
expressions for the U-matrices. For R <^ L, H and 
arbitrary H/L the result for the force can be written as 
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The fj can be computed exactly. We have obtained them 
up to j ~ 11 and the first three are (with s = \/l + 4/1^) 
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The coefficient fr of vanishes since a multipole of or- 
der I contributes to the T-matrix at order so that 
beyond the two-dipole term ~ the next term comes 
from a dipole (1=1) and a quadrupole {I = 2), yielding 
/g. For > L one has fe{h) = -1001/16 + 3/{Ah^) + 
Oih-^), fsih) = -71523/160 + 39/(80/i6) + ©(/i-S) so 
that the wall induces weak repulsive corrections. For 
H <t: L, fe{h) = -791/8 + 6741/iV8 + ©(/i^), fs{h) = 
-60939/80 + 582879/iV80 + 0{h^) so that the force am- 
plitude decreases when the spheres are moved a small 
distance away from the wall. This proves the existence 
of a minimum in the force amplitude as function of H/R 
for fixed, sufficiently small R/L. We note that all fj{h) 
are finite for h ^ oo but some diverge for h 0, e.g., 
/g ~ fii ^ , making them important for small H. 

To obtain the interaction at smaller separations or 
larger radius, we have computed the energy £"00 and force 
F = —d£oo/dL between the spheres numerically. For the 
energy, we have computed the last determinant of Eq. ([3| 
and the integral over k of Eq. ([T|) numerically. The force 
is obtained by polynomial interpolation of the data for 
the energy. The matrices are truncated at a sufficiently 
large number of partial waves (with a maximum trunca- 
tion order Imax = 17 for the smallest separation) so that 
the relative accuracy of the values for £00 is ~ 10"'^. The 
data for H/R = 1 are obtained by extrapolation in Imax- 
The results are shown in Figs. |3] |4j In order to show the 
effect of the wall, the figures display the energy and force 
normalized to the results for two spheres without a wall. 
Fig. [3] shows the energy and force as function of the (in- 
verse) separation between the spheres for different fixed 
wall distances. Energy and force show an increasing rela- 
tive enhancement due to the wall with increasing L, with 
the maximal enhancement for small H . For sufficiently 
large H/R, the energy and force ratios are non-monotonic 
in L and can be slightly smaller than without the wall. 
Fig. |4]shows the force as function of the wall distance for 
fixed L. When the spheres approach the wall, the force 
first decreases slightly if R/L < 0.3 and then increases 
strongly under a further reduction of H. For R/L > 0.3 
the force monotonically increases towards the wall. This 
agrees with the prediction of the large distance expan- 
sion. The expansion of Eq. ( 11 1 with j = 10 terms is also 
shown in Fig. |4]for R/L < 0.2. Its validity is limited to 
large L/R and not too small H/R; it fails completely for 
R/L > 0.2 and hence is not shown in this range. 

Our results for atoms are potentially relevant to the 
interaction between trapped Bose-Einstein condensates 
and a surface at close surface separations. The results 
for macroscopic spheres could be important for the de- 
sign of nano-mechanical devices where small components 
operate in close vicinity to metallic boundaries. Gen- 
erally, the wall-induced enhancement of the interaction 
can make the experimental observation of Casimir forces 
between small particles more feasible. The reported de- 
pendence on the anisotropy of polarizabilities suggest in- 



teresting effects for objects of non-spherical shape. 
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FIG. 3: Numerical results for the potential energy (dashed 
curves) and force (solid curves) between two spheres as func- 
tion of R/L for different sidewall separations H/R. Both force 
and energy are normalized to their values in the absence of 
the sidewall. Inset: Magnification of behavior for small R/L. 
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FIG. 4: Numerical results for the force (dots) between two 
spheres as function of the sidewall separation H/R for dif- 
ferent sphere separations R/L. Shown are also the analytical 
results of Eq. (Ill, including terms up to j = 10 for R/L < 0.2 
(solid curves). Inset: Magnification of the non-monotonicity. 
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